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Abstract 

Kinematic fitting, where constraints such as energy and momentum conservation are 
imposed on measured four-vectors of jets and leptons, is an important tool to improve the 
resolution in high-energy physics experiments. At future e + e~ colliders, photon radiation 
parallel to the beam carrying away large amounts of energy and momentum will become 
a challenge for kinematic fitting. A photon with longitudinal momentum p zri (rj) is in- 
troduced, which is parametrized such that rj follows a normal distribution. In the fit, rj is 
treated as having a measured value of zero, which corresponds to p Zj7 = 0. As a result, fits 
with constraints on energy and momentum conservation converge well even in the presence 
of a highly energetic photon, while the resolution of fits without such a photon is retained. 
A fully simulated and reconstructed e + e~ — > qqqq event sample at ^/s = 500 GeV is used 
to investigate the performance of this method under realistic conditions, as expected at the 
International Linear Collider. 



To be submitted to Nucl. Instrum. Meth. A 



1 Introduction 



Radiation of photons at angles so small that they escape along the beam pipe is usually not 
taken into account in kinematic fits. At previous e + e~ colliders such as LEP, the losses due 
to photon radiation were acceptable [1]. At future facilities such as the International Linear 
Collider (ILC) or the Compact Linear Collider (CLIC) 0, photon radiation will be much 
stronger due to higher center-of-mass energies and stronger focussing of the beams, which 
makes it desirable to model photon radiation in kinematic fits. 

Kinematic fitting is a well-established tool to improve jet energy and invariant mass resolu- 
tions. A number of four- vectors representing the final state particles is fitted under constraints 
such as energy and momentum conservation. The four-vectors are parametrized by suitably 
chosen variables such that the measured values follow an approximately Gaussian distribution 
around the true values. A \ 2 that quantifies the deviation between measured and fitted parame- 
ters is minimized under the condition that the imposed constraints are fulfilled flU. 

The improvement in resolution emerges from the redundant information contained in the 
measured values in the presence of constraints. Unmeasured parameters reduce the redundancy, 
since one constraint is used up for each unmeasured parameter to determine its value. The 
redundancy is quantified by the number of degrees of freedom, which is given by the number 
of constraints minus the number of unmeasured parameters. 

The two main effects that cause the emission of photons collinear with the incoming beams 
so that they escape the main detector are initial state radiation (ISR) and beamstrahlung. ISR 
is a higher-order QED effect, at which real photons are emitted before the actual interaction. 
Beamstrahlung is caused by the electrical fields of the bunches colliding with each other: elec- 
trons in the one bunch are deflected by the field of the other bunch and thus emit bremsstrahlung 
photons. 

ISR is characterized by an energy spectrum that follows a power law with an exponent of 
roughly —0.9 0. Thus the vast majority of events have at most one ISR photon with an energy 
above a few GeV, which is the accuracy to which the total energy and longitudinal momentum 
of fully hadronic events can be measured by a typical detector envisioned for the linear collider. 
This photon can, however, carry substantial energy of tens of GeV. Beamstrahlung on the other 
hand has an energy spectrum with an exponent of —2/3, but with an additional exponential 
suppression of high energy photons 10. The mean number of beamstrahlung photons emitted 
prior to the interaction can be of order one or even larger, depending on the beam parameters. 

This paper presents a novel method to take the energy and longitudinal momentum of photon 
radiation into account in kinematic fits. A priori information about the momentum spectrum of 
photon radiation is used to treat the photon's momentum as a measured parameter in the fit. 
As a test case, the production of W + W~/Z°Z° pairs decaying to light quark jets at the ILC is 
considered, with fully simulated Monte Carlo events as reconstructed by the International Large 
Detector (ILD) Q. A more detailed description of the method and its application tests can be 
found in 0. 

The main focus of this method is an improved treatment of the effects of ISR, because ISR 
is the main source of highly energetic photons. Therefore, only a single photon is included 
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in the kinematic fit, with an energy spectrum given by a power law, as expected for ISR. A 
similar method with the inclusion of two photons in the fit and an energy spectrum describing 
the combined effects of ISR and beamstrahlung is the subject of ongoing work and beyond the 
scope of the current publication. However, the method presented here leads to a significant 
improvement also in the presence of beamstrahlung, as shown in section[3] 



2 Representation of the photon 

Since photons from ISR and beamstrahlung escaping the detector have to a good approximation 
zero transverse momentum with respect to the beam direction, they affect mainly the conser- 
vation of (detected) energy E and longitudinal momentum p z . The simplest method to cope 
with highly energetic photons in a constrained kinematic fit is therefore to drop the energy and 
longitudinal momentum conservation constraints, thus losing two degress of freedom. 

A somewhat better solution is to introduce a fit object representing the undetected photon 
with one free, unmeasured parameter, namely its longitudinal momentum p z>y , and set p xn = 
Py i7 = and thus E y = \p Zj7 \. This allows the energy and p z constraints to be recovered, at the 
price of one unmeasured parameter, so that one degree of freedom is regained. 

However, this approach neglects the information about the momentum spectrum of the pho- 
tons. Here this information is used so that the photon is treated as a particle with a measured 
momentum of zero and an uncertainty derived from its known momentum spectrum. 

2.1 Parametrization of the Photon Energy 

In a kinematic x 2 fit the measured four- vector components of a particle or jet are parametrized 
with parameters r/i (e.g., E, 9, <p) such that the difference ?7i, m eas — ?7i,tmc between the measured 
?7i,mcas and the true value ?7i itmc follows a Gaussian distribution with zero mean and standard 
deviation 5rji (for reasons of notational simplicity we limit the discussion to the case where the 
parameters r\\ are uncorrected). Then, 



is, apart from a constant, identical to the negative logarithm of the likelihood to obtain the 
measured values, given the values r\\. 



Thus, the \ 2 fit seeks the best estimate r\\ for the true parameter values by maximizing the 
likelihood to get the observed parameter values r/ mca s under the condition that the imposed 
constraints are fulfilled, which are expressed by a number of constraint functions gy (tji) = 0. 
No assumption is made, or is necessary, about the distribution of the true parameter values 




(1) 



X 2 = - InV (?7i, m eas|?7i) + const. 



(2) 
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However, if an ensemble of events is considered where the distribution of a parameter ?7 truo 
is known to be Gaussian with zero mean, then for this ensemble the choice 77 mcas = also leads 
to a Gaussian distribution of ?7 meas — r/tmo and for such an ensemble it is justified to estimate 
?7true by means of a x 2 fit- 

In the case of photon radiation, the distribution of the unmeasured momentum p zr/ is known, 
though definitely non-Gaussian. Thus we seek a parametrization of the photon's momentum 
Pz,7 = Pz,7 (v) sucn that the true value of r] follows a Gaussian distribution with mean zero and 
unit standard deviation 8rj — 1. Then the photon will be treated as if it had a measured value of 
f]raeas = 0. The photon will then be added to the list of fit objects in the kinematic fit, thereby 
introducing an additional contribution to the overall \ 2 °f f] 2 /^ 7 ] 2 = V 2 - By this procedure, the 
a priori knowledge of the photon's energy spectrum (in particular the fact that it is negligibly 
small in most cases) is used, and all energy and momentum constraints can be applied. 

The probability density function V (y) for the energy fraction y = E 7 /E heaxa carried by 
initial state radiation is well approximated by [0 



with the exponent given by 



P(y)=Py f, -\ (3) 



7T \ ?7ig 



which corresponds to (3 = 0.1235 for y/s = 500 GeV. 

Considering that an ISR photon can be emitted by either beam leads to 

P 



z,7 J 



2E n 



EL 



(5) 



where E max < -Ebcam is the maximum possible photon energy. As a consequence, the quantity 
z given by 

z = sign(p z , 7 ) (6) 



max 



is uniformely distributed between —1 and I, and hence 

rj = v^-erf" 1 ^), (7) 
follows a Gaussian distribution with zero mean and unit standard deviation. Here, erf -1 (z) 

X 

denotes the inverse of the error function given by erf (x) — j e~* dt. 

n o 

Conversely the expressions for z and p z j7 as a function of the parameter 77 read 

z( V ) = erf(77/\/2) (8) 
Pz^iv) = sign(z) E m ^ x \z\P (9) 

= sign(?7) £ max ed(\r]\/V2) 



* • (10) 
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2.2 Properties of the Parametrization 



0.1235. The function has four 
= 0, reflecting the fact that the 
77 1 > 0.7 significant momenta 

Around 77 = the value of p zr/ does not change and thus cannot influence the global x 2 
of the kinematic fit. Therefore, the penalty term rf leads to a local minimum of the x 2 at this 
value of 77. This will also be the global minimum if the measured four-momenta of the final 
state particles are compatible with no missing momentum from ISR. In this case, a fit with a 
photon fit object has exactly the same result as a fit without a photon. 

Due to this local minimum, any minimization method based on derivatives will always yield 
T] = if this value is used as starting value in the minimization. Therefore, to find the global 
minimum, in addition a different starting value must be tried, for instance 77 (— p z . m iss) calculated 
from the missing p z of the event. 

For 1 < |7;| <2.5, the curve rises steeply, so that large values of missing energy and mo- 
mentum from ISR and beamstrahlung can be accomodated by the kinematic fit at a moderate 
penalty rf. Thus, for photon momenta that are large compared to the detector resolution, the 
kinematic fit should find a global minimum of the x 2 close to the true photon momentum, with 
a negligible bias towards low photon energies, despite the fact that the "measured" value of rj 
and thus of p z 7 is set to zero. 

Above \r)\ ~ 2.5, corresponding to |p Z)7 |/-E m ax ~ 0.9, the curve flattens again, so that ex- 
tremely large photon momenta are suppressed due to the fact that they are increasingly unlikely. 
This region is, however, of little interest in realistic analyses. 



3 Performance Tests 

The method described above is applied to the process e + e~ — > W + W~ — > 4jets events. The 
fraction of successful fits, the width and the shift of the reconstructed W ± mass peak are used 
to compare the performance of the various kinematic fit variants. 

3.1 Data Set 

The analysis sample was generated using the matrix element generator WHIZARD [0, which 
takes into account all Feynman diagrams leading to a given final state, including interference 
terms. Here, the process e + e" — > uddu is chosen, because it contains no heavy quarks in the 
final state, so that the jet energy and angle measurement is not compromised by the presence of 
neutrinos from semileptonic decays. 

In addition to the dominant W ± pair production, also Z° pair production contributes to the 
formation of this final state, and at a even smaller level single boson production with subsequent 



Fig. \T\ shows a graph of p Z)7 (77) for E m3iX = 225 GeV and /3 = 
distinct bends around \rj\ ~ 1 and \rj\ ~ 2.5. It is flat around 77 
majority of ISR photons have negligible momentum; only for 
above 1 GeV are predicted. 
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radiation off final state quarks. Due to the inclusion of interference effects, it is conceptually 
not possible to identify events from these additional processes and remove them. The centre- 
of-mass energy is y/s = 500 GeV, and the W ± mass was set to m^ n = 80.419 GeV. 

The initial state radiation is also simulated by WHIZARD, and one ISR photon per incom- 
ing lepton is stored in the event record. In contrast, the energy spread of the incoming beams 
and beamstrahlung is taken into account in the event generation by a corresponding variation 
of the momenta of the incoming leptons, using a beamstrahlung spectrum that was simulated 
with GUINEA-PIG ifTOI . For this calculation the nominal beam parameter set of the ILC was 
assumed, ifTTTl . in particular an energy spread of 0.14 % and 0.07 % for the electron and positron 
beams, respectively, a beamstrahlung parameter T ave = 0.047, a mean energy loss by beam- 
strahlung of 5bs = 0.023, and a vertical disruption parameter of D y = 19.1. 

Thus, the momenta of ISR photons are directly accessible in the event record, while the 
combined effect of beamstrahlung and beam energy spread has to be deduced from the total 
four-momentum of all final state particles. 

A full simulation of the ILD detector Q is performed by the GEANT based simulation pro- 
gram MOKKA lfi~2l . In the event reconstruction, which is implemented as part of the software 
package MarlinReco lfi~3~1 . the tracks are matched to the calorimeter clusters by the Pandora par- 
ticle flow algorithm [fT4l and the resulting reconstructed particles are forced into four jets by the 
Durham algorithm lfl~5l . Each of the four jets has to have a minimum energy of E^ ct > 5 GeV 
and a polar angle that fulfills | cos #j et | < 0.9. 

The jet momentum four- vectors are parametrized in terms of energy E- iCt , polar angle 9 and 
azimuthal angle cf), with resolutions |[T6ll : 

5Ej Ct /E ict = 32.24 %/^/E~ t + 1.242 -lO" 4 ^- 1.446% (11) 

60 = 0.03925/ v /^ + 0.3373/E jct (12) 

50 = 0.05873/v/!~^ + 0.3207/£ jet . (13) 

Since the method presented here concerns radiation escaping the main detector, a subsample 
of events is selected such that at generator level only negligible energy from ISR is present in the 
detector acceptance. Therefore, events are rejected with ISR photons of energy E& cn > 5 GeV 
and polar angles 0.29° < 6^ cn < 179.71°, which corresponds to the acceptance of the beampipe 
calorimeter (BeamCAL). No cut is applied on the energy or direction of the beamstrahlung. 

In order to investigate the influence of ISR and beamstrahlung on the performance of the 
kinematic fit, the event sample is divided into three subsamples according to the total energy 
Ei SR of the ISR photons: 

• A subsample with small -Eisr < 5 GeV, where the ISR is expected to have a small effect 
and thus a kinematic fit is expected to perform well without an additional photon, is used 
to evaluate whether the addition of such a photon leads to a loss of resolution. 

• A subsample with moderate ISR energy (5 < -Eisr < 30 GeV) is used to evaluate a 
possible bias of the kinematic fit, and whether the addition of a photon removes this bias 
and increases the fraction of good fits. 

• A subsample with large Ej SR > 30 GeV serves to quantify how well the fit with photon 
works, and whether it has any advantage over a fit where the energy and longitudinal 
momentum constraints are dropped completely. 



5 



3.2 Evaluation Method 



In order to investigate the performance of the proposed method, kinematic fits are applied to 
the four jets in the events from the test sample, comparing the event hypotheses "4 jets" (4j) 
and "4 jets + 1 photon" (4j + 7). Both event hypotheses are fitted with five constraints (5C- 
fit): conservation of energy, conservation of the three momentum components and equal di-jet 
masses. In addition, the events are fitted also using only the three constraints (3C-fit) that are 
not affected by the presence of photon radiation, i.e. conservation of the transverse momentum 
components and the equal mass constraint. 

As discussed in Sect. 12. 1 1 both values p z j7 = and p zn = — p Zjm - lss are considered as starting 
values for the photon momentum in the kinematic fit, and the result with the better x 2 is chosen. 
Values of (3 = 0.1235 and Eisr = 225 GeV, which is the maximal photon energy that allows 
W ± pair production, are used in the photon parametrization Eq. (flOl) . Fig.[2]shows the quantities 
z and rj of Eqs. © and ©, calculated from p z i7 of the most energetic ISR photon in the event. 
It can be seen that indeed z is distributed uniformly and rj follows a Gaussian distribution. 

The sample used for the performance tests includes the effects of both ISR and beam- 
strahlung. Since the photon parametrization used here has been derived from the ISR momen- 
tum spectrum, tests are first performed that exclude the effect of beamstrahlung. The beam- 
strahlung is artificially "turned off" by using the total generated energy and momentum of the 
final state particles, including the ISR photons, in the energy and momentum constraints, rather 
than the nominal values of Px, y ,z = 0, E = = 500 GeV. Alternatively, the constraints 
are set to these nominal values, so that the combined effects of ISR and beamstrahlung can be 
studied. 

An important indicator for the performance of the various kinematic fits considered is the 
fraction of good fits, which are defined as those having a fit probability p > 0.001. 

Due to the intrinsic widths of the W ± and Z° bosons, which are not negligible compared to 
the detector resolution at an ILD-type detector, the equal-mass constraint is only approximately 
fulfilled by the four- vectors on generator level and therefore reduces the fraction of good fits. 

The equal-mass constraint is applied mainly in order to choose the correct jet pairing. The 
fit is performed for all three possible jet pairings, and the pairing that results in the best x 2 value 
is chosen as the correct one under the assumption that the jets stem from either a W + W~ or 
a Z°Z° pair. Because of the W ± and Z° width, the equal-mass constraint leads to an aver- 
age x 2 contribution that is significantly higher than the value expected for the addition of one 
constraint, i.e. one degree of freedom. As a consequence, only 55 % of the 3C-fits have a fit 
probability p > 0.001. A more elaborate treatment of the equal-mass constraint that would 
increase this fraction is, however, beyond the scope of the present analysis. 

Fig.[3]shows the the invariant di-jet masses before and after the kinematic fit for the complete 
sample, including ISR and beamstrahlung. A clear peak at the W ± mass is observed, while the 
much smaller Z° mass peak appears only as an enhancement on the right side of the W ± mass 
peak. 

Imposing an equal-mass constraint leads to an implicit averaging of the two di-jet masses 
in each event. Therefore, the average di-jet mass before the fit is compared to the dijet masses 
after the various kinematic fits in Fig.[3l 
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For a quantitative comparison of the different kinematic fits, each mass distribution is fitted 
with an analytic function. The W ± mass peak is expected to follow a relativistic Breit-Wigner 
distribution, folded with a Gaussian distribution that reflects the detector resolution. Here, the 
mass peaks are fitted with a Voigt function V a ^{x) ifTTl . which is the convolution of a non- 
relativistic Breit-Wigner (Cauchy) distribution of width Y and a Gaussian distribution with an 
RMS of a. Thus the following function is fitted to the histograms in the range 75 < m < 
95 GeV: 

f(m) = N - ((1 - f z ) ■ V a>rw (m - m w ) + fz ■ V a> r z {m - m z )) (14) 

The values for the Z° mass mz = 91.19 GeV and the decay widths^ T\y = 2.14 GeV and 
Tz = 2.50 GeV are fixed to their literature values [Tl8l . The same Gaussian width a, reflecting 
the detector resolution, is used for the W ± and Z° mass peaks. N corresponds to the number of 
histogram entries and f z to the fraction of Z-pair events. However, because the symmetric Voigt 
function does not describe the asymmetry of a relativistic Breit-Wigner distribution correctly, 
fz is not an accurate estimate of the fraction of Z° events in the sample. 

The parameters of interest are the Gaussian width a and the difference Am w = m w — m^ n 
between the fitted W ± mass m w and the input W ± mass mf^ n . 

If large amounts of energy are missing, the fitted jet energies have to be larger than the mea- 
sured ones to fulfill energy conservation. Consequently, di-jet masses are shifted to higher val- 
ues and thus a larger Am w is obtained. Due to the imperfections of the lineshape fit, a nonzero 
value of Am w is to be expected, for which a correction would be applied in a real analysis. 
However, if this mass shift depends on the amount of energy from ISR and beamstrahlung, it 
leads to a broadening of the signal and thus a loss of resolution; in addition, systematic un- 
certainties arise from the description of the ISR and in particular the beamstrahlung energy 
spectrum. Therefore, a mass shift that is independent of the amount of energy lost to ISR and 
beamstrahlung is desirable. 

3.3 Results 

Tab. \T\ summarizes the results of our tests. It lists the fraction of good fits, the mass shift and the 
width of the Gaussian part of the Voigt function for the complete sample, as well as the three 
subsamples with different amounts of missing energy due to ISR photons. The results are given 
for the average of the di-jet masses before a kinematic fit, using the 3C jet pairing, as well as 
the di-jet mass after applying a 3C fit or a 5C fit without or with an ISR photon. The results are 
reported for the case where the effect from beamstrahlung has been excluded by adjusting the 
energy and momentum constraints (cf. Sect. 13.21) , and for the realistic case where effects from 
ISR and beamstrahlung are fully taken into account. 

Results with ISR only 

A comparison of the fit results demonstrates the gain in resolution achieved by kinematic fitting: 
The Gaussian a, which corresponds to the di-jet mass resolution, is a = 2.1 GeV for the 

1 The average of two independent random numbers distributed according to a Breit-Wigner of width V follows 
a Breit-Wigner distribution of the same width T. 
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Sub sample 
(Fraction) 


Constraints, 
Hypothesis 


ISR only 
Good Amw cw 
fits [%] [GeV] [GeV] 


Full Photon Spectrum 
Good Amw c"w 

t-% r n r . — -r T i r s~--i -r -r~\ 

fits [%] [GeV] [GeV] 


All events 

(100%) 


3C,4j 
5C,4j 
5C,4j + 7 


55% +0.78 2.05 
55% +0.82 2.06 
42% +0.67 1.21 
54% +0.53 1.25 


55% +0.78 2.05 
55% +0.82 2.06 
31% +0.91 1.30 
52% +0.75 1.35 


E 1SR < 5 GeV 
(75 %) 


3C,4j 
5C,4j 
5C, 4j + 7 


56% +0.80 2.04 
56% +0.85 2.06 
53% +0.63 1.19 
55% +0.49 1.24 


56 % +0.80 2.04 
56% +0.85 2.06 
40% +0.86 1.27 
54% +0.69 1.31 


5 < E ISR < 30 GeV 
(11%) 


3C,4j 
5C,4j 
5C, 4j + 7 


54% +0.79 2.07 
54% +0.84 2.08 
15% +1.68 1.25 
53% +0.71 1.27 


54% +0.79 2.07 
54% +0.84 2.08 
12% +2.19 1.29 
50% +1.07 1.51 


£isr > 30 GeV 
(13%) 


3C,4j 
5C,4j 
5C, 4j + 7 


53% +0.59 1.99 
53% +0.66 1.99 

0% 

47% +0.64 1.21 


53% +0.59 1.99 
53% +0.66 1.99 
0% 

42% +0.91 1.38 



Table 1 : Results of kinematic fits under various conditions. "ISR only" refers to the case where 
the effect of beamstrahlung and beam energy spread is removed from the fit as explained in the 
text, while "Full Photon Spectrum" includes these effects. For each fit variation, the fraction of 
good fits with fit probability p > 0.001, the difference Amw between the fitted and generated 
W mass of m^ n = 80.419 GeV, and the width of the Gaussian part of the Voigt function is 
given. The rows refer to the results from averaging the measured di-jet masses without a fit 
for events where the 3C fit converges, the 3C fit with only transverse momentum and equal- 
mass constraint, the 5C fit under a four jet hypothesis with longitudinal momentum and energy 
constraints in addition, and the 5C fit with an additional ISR photon fit object. The subsamples 
are distinguished by the total energy Ei S r of ISR photons, excluding beamstrahlung. 

average of the two di-jet masses without a kinematic fit and improves to a = 1.3 GeV if a 
kinematic fit with five constraints is used. A fit with only three constraints does not improve the 
resolution compared to the simple averaging of the unfitted di-jet masses. 

The fit with five constraints and no ISR photon cannot be applied to the subsample with 
Eisk > 30 GeV, because fit probabilities above the cut of p = 0.001 are essentially never 
achieved due to the missing energy and momentum that are are too large to be accomodated 
by the experimental resolution of a few GeV. Therefore this subsample, which contains 13 % 
of all events, cannot be used for an analysis. The 5C fit with an ISR photon, on the other 
hand, achieves almost the same performance for the two subsamples with Eisr > 30 GeV and 
-E'isr < 5 GeV in terms of the fraction of good fits (47 % vs. 55 %) as well as in resolution (a = 
1.21 GeV vs. 1.24 GeV) with only a small additional bias in the W mass (Am w = 0.64 GeV 
vs. 0.49 GeV). 

The sample with moderate ISR energy 5 < Eisr < 30 GeV, which comprises 1 1 % of the 
events, demonstrates that the 5C fit without the inclusion of an ISR photon tends to develop 
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a mass bias. This is because the energy carried away by the photon is falsely attributed to 
the final state jets, which increases their energy and thus the invariant mass: The mass bias 
increases from Am w = +0.63 GeV to +1.68 GeV. At the same time, only 15 % of the events 
yield a good 5C fit under the Aj hypothesis. In contrast, the Aj + 7 hypothesis shows the same 
performance in terms of fraction of good fit, mass shift and resolution as for the sample with 
small missing energy. 

The fact that for all fit hypotheses only about half of the events have reasonable fit proba- 
bilities p > 0.001 can be mostly attributed to the equal-mass constraint: The resolution for the 
difference of the di-jet masses is approxiately 4.1 GeV (twice the resolution for the di-jet mass 
average for the unfitted jets), which is of similar size as the broadening^ of 4.3 GeV due to the 
intrinsic W width. This indicates that in a real analysis the naive equal-mass constraint has to 
be modified to take the natural W width into account. Other factors that reduce the fraction 
of successful fits are events from processes other than W/Z boson pair production and the fact 
that the jet error parametrization employed in this analysis does not include the effects of parton 
showering. 

Results with ISR and beamstrahlung 

The right-hand side of Tab. \T\ shows the results for the case where the effect of both, ISR and 
beamstrahlung, is considered. Because the three subsamples are defined on the basis of the ISR 
energy only, the same amount of beamstrahlung is present in each of them. A comparison with 
the case where only the effect from ISR is considered, demonstrates that the photon momentum 
parametrization Eq. ( flOl) derived from the ISR momentum spectrum also works quite well in 
the presence of beamstrahlung, at least at the level of beamstrahlung that is expected for the 
nominal ILC parameter set. 

Since beamstrahlung in the Monte Carlo simulation used for this analysis is simulated solely 
through a variation of the energy of the incoming leptons, no transverse momentum is carried by 
the beamstrahlung. Therefore the results for the 3C fit and the di-jet masses calculated without 
a kinematic fit do not change when beamstrahlung effects are considered. 

The performance of the 5C fit under the Aj hypothesis is significantly reduced when beam- 
strahlung effects are considered due to the larger amount of missing energy. Overall, the frac- 
tion of good fits goes down from 42 % to 31 %. For the subsample with less than 5 GeV of 
ISR energy it is reduced from 53 % to 40 %. At the same time, the W ± mass shift increases by 
approximately 0.2 GeV for the whole sample. For the subsample with medium Ei SK , however, 
the mass shift increases from +1.68 GeV to +2.19 GeV. 

On the other hand, with the Aj + 7 hypothesis, the 5C fit performance is much less affected 
by beamstrahlung effects: The fraction of good fits stays almost constant, and the a of the 
Gaussian width of the mass peak increases only moderately, from 1.25 GeV to 1.35 GeV for 
the complete sample. The mass shift increases by approximately 0.2 GeV for the full sample, 
which is similar to the Aj hypothesis. However, for the subsample with 5 < Ei SR < 30 GeV 
the mass shift is significantly reduced from +2.19 to +1.07 GeV by the inclusion of the photon 

2 The difference of two Breit-Wigner distributed random numbers follows itself a Breit-Wigner function with 
a width that is the sum of the two individual widths. 
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in the fit. The increase of the mass shift with respect to the ISR only case indicates that the 
Aj + 7 hypothesis cannot fully accomodate beamstrahlung effects, because typically both beam 
particles radiate off significant energy. This may necessitate the inclusion of a second photon in 
the fit. 

As a final check, Fig. |4] shows the fitted longitudinal momentum p z r/ of the photon versus 
the generated pf c ° of the most energetic ISR+beamstrahlung photon pair in the event, where the 
momenta of the ISR and beamstrahlung photons with either positive or negative p z are added. It 
can be seen that the fitted photon momentum p za corresponds quite well to the true momentum, 
without any visible bias. In particular, the fact that the photon is treated as having a measured 
p Z)7 = does not lead to a large bias towards small values of p zr/ . This is explained by the fact 
that the function p zr/ (77) of Eq. (flOl) rises very rapidly, as discussed in Sect. 12.11 

The right side of Fig |4] shows the difference Ap Zj7 = sign(p Zi7 ) • (p Zi7 — pf™). The mean 
(Ap Z 7 ) = —0.32 GeV is small, and negative, showing that the reconstructed |p z 7 | is slightly 
smaller on average than the generated one, as expected, but that this bias is indeed quite small. 
The resolution for p zrf is found to be 3.25 GeV. 



4 Summary and Conclusions 

In this paper a method is proposed to take the effect of ISR into account in kinematic fits by 
introducing a photon that is treated as if its measured momentum were zero. The longitudinal 
momentum p z l is expressed as a function p z l (77) of the parameter 77 such that the true value of 
7] follows a normal distribution with zero mean and unit standard deviation. 

The performance of this method is evaluated using a sample of e + e _ — > uddu events, which 
is dominated by W + W~ pair production, at y/s = 500 GeV. The sample includes the effects 
from ISR and beamstrahlung. It is fully simulated and reconstructed, using the simulation 
for the ILD detector at the ILC. A 5C kinematic fit with energy and momentum conservation 
constraints and an equal-mass constraint is applied, and the results for the fit hypothesis with 
four jets and a photon are compared to three alternatives: a 5C fit with a conventional four jet 
hypothesis, a 3C fit where the energy and longitudinal momentum constraints are dropped, and 
the results obtained without a kinematic fit. 

The 5C fit with the new Aj + 7 hypothesis performs as well as a 5C fit with a Aj hypothesis 
in terms of resolution, while a 3C is significantly worse and does not yield any improvement 
over a mass reconstruction without any kinematic fit. 

For events with significant energy from ISR photons (5 < -Eisr < 30 GeV), the fraction of 
good fits with a fit probability p > 0.001 drops from 40 % to 12 % for a 5C fit without a photon, 
and goes to zero for Eisr > 30 GeV. In addition, as the missing energy is distributed to the 
jets by such a fit, a shift of the reconstructed di-jet masses towards larger values is observed. 

Both problems are solved by the new Aj + 7 hypothesis: even for large values of Ei SR > 
30 GeV, the fraction of good fits and the di-jet mass resolution are similar to the values obtained 
at Eisr < 5 GeV, while the mass shift remains small. 
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In short, under the Aj + 7 hypothesis, a 5C fit achieves the same resolution as with a con- 
ventional Aj fit hypothesis, but independent of the amount of ISR energy, without developing a 
mass bias, and with a similar fraction of good fits as a 3C fit. 

Although the parametrization p z a (77) was developed using the momentum spectrum of ISR 
photons, the method also performs well in the presence of beamstrahlung, at least at the mod- 
erate level expected for the nominal parameter set of the ILC. 

In a future development the parametrization could be adapted to include beamstrahlung 
effects. This may be necessary in scenarios with enhanced beamstrahlung, such as the "low 
power" parameter set proposed for the ILC, or at CLIC. We expect that under such conditions 
the addition of a second photon in the fit would become necessary in order to take into account 
the energy loss suffered by both beam particles. 
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Figure 2: The parameters z (a) of Eq. © and 77 (b) of Eq. ©, calculated from the longitudinal 
momentum p z j7 of the most energetic ISR photon in the Monte Carlo sample described in the 
text, using E max = 225 GeV and (3 = 0.1235. z is expected to be uniformely distributed in 
— 1 < z < 1, and 77 should follow a Gaussian distribution with zero mean and unit standard 
deviation, which is shown for comparison in the plot. 
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Figure 3: Invariant di-jet masses for the Monte Carlo sample described in the text: a) the 
average of the two di-jet masses calculated from the measured four-vectors, using the jet pairing 
from the 3C fit sample; b) rriy s after application of the 3C fit; c) my for the 5C fit under a 4j 
hypothesis; d) for the 5C fit under a 4 j + 7 hypothesis. 




C[ QeV ] Ap zY [GeV] 



Figure 4: Fitted photon momentum p za plotted against the true momentum pf? of the most 
energetic ISR+beamstrahlung photon combination in the event (a), and the difference Ap z i7 = 

sign(p Zi7 ) ■ (p Zi7 - pf^ 1 ) (b). 
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